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The superfluid density p s in underdoped (T c ~ 23 K) , optimally doped (T c ~ 35 K) and overdoped 
(T c ~ 29 K) single crystalline (BiPb)2(SrLa)2Cu06+« samples was studied by means of muon-spin 
rotation (^iSR). By combining the ^tSR data with the results of ARPES measurements on similar 
samples [Nature 457, 296 (2009)] good self-consistent agreement is obtained between two techniques 
concerning the temperature and the doping evolution of p s . 
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The superfluid density p s , being proportional to the 
density of the supercarriers, is one of the important char- 
acteristic of the superconducting condensate. Consider- 
ing that superconductivity is characterized by the phase 
coherence of electrons forming the pairs (two-particle 
process), it was previously believed that techniques 
which probe the single-particle excitations of the con- 
densate, such as photoemission and single-electron tun- 
neling, could not directly provide the information on p s . 
It was quite unexpected, therefore, when angle-resolved 
photoemission (ARPES) studies of the cuprate high- 
temperature superconductor (HTS) Bi 2 Sr 2 CaCu 2 08+(5 
revealed that a sharp peak, formed below the supercon- 
ducting transition temperature T c near the Brillouin zone 
boundary, contains information not only on the pairing 
strength (the superconducting energy gap A) but also on 
the phase coherence [j], 0] . The peak intensity shows a 
clear resemblance to the behavior exhibited by Ps and 
scales linearly with T c in the underdoped regime 

E li- 
lt should be mentioned that previous comparisons be- 
tween p s and the coherence peak (CP) were made for CP 
measured near the Brillouin zone boundary 0, Q ■ The 
superfluid density, in its turn, is an "angular integrated" 
quantity accumulating information over the full Fermi 
surface. More importantly, the carriers near the antin- 
odes could be affected by the pscudogap, which within 
"two-gap" scenario is supposed to be unrelated to the 
This could lead to an additional decrease 
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pairing 

of the CP intensity in the antinodal region ll( . Conse- 
quently the CP measurements over the full Brillouin zone 
are needed in order to compare ARPES data with the su- 
perfluid density studied independently in e.g., muon-spin 
rotation (/iSR) or microwave experiments. 

In this paper we report on the results of pSK stud- 
ies of the superfluid density and its comparison with 
the previously reported ARPES data [ll| for underdoped 
(T c ~ 23 K), optimally doped (T c ~ 35 K) and overdoped 



(T c ~ 29 K) single crystalline (BiPb) 2 (SrLa) 2 Cu0 6+( 5 
samples. It was found that p s (T) could be well described 
with a superconducting gap of d— wave symmetry by 
introducing the coherence quasiparticle weight function 
QW measured by means of ARPES. The T dependence 
of the superconducting energy gap A follows the BCS 
prediction, and p s (T = 0) scales with the CP intensity 
integrated over the whole Fermi surface. 

Details on the sample preparation of 
(BiPb) 2 (SrLa) 2 Cu0 6 +A (Bi2201) single crystals can 
be found elsewhere jllMl3| . The samples are labelled 
by their superconducting transition temperature T c 
with the prefix UD for underdoped, OP for optimally 
doped, and OD for overdoped as: UD23K, OP35K, and 
OD29K. Zero-field (ZF) and transverse-field (TF) ^tSR 
experiments were carried out at the 7rM3 beam line 
(Paul Scherrer Institute, Switzerland). The samples 
were cooled from above T c to 1.6 K at H = during 
ZF-/iSR experiments and in a series of fields ranging 
from 5 mT to 0.64 T in TF-/^SR experiments. In TF 
studies the magnetic field was applied parallel to the c 
axis and transverse to the muon-spin polarization. The 
typical counting statistics were ~ 15 — 18 million muon 
detections per data point. The TF-uSR data for OP35K 
are partially published in Refs. [3 UU 

The representative ZF and TF-/^SR time-spectra mea- 
sured above and below T c are presented in the left panels 
of Figs.[T]a, b, and c. In ZF all three samples show a slow, 
temperature independent decay thus implying that the 
T dependent relaxation observed in TF experiments at 
T < T c should be attributed to the inhomogencous field 
distribution in a superconductor in the vortex state. The 
ZF data for OP35K and OD29K are well described by 
the Gaussian decay function A ZF (t) = A exp(— a^t 2 /2) 
which is caused by the dipolar field arising from the nu- 
clear magnetic moments (Aq is the initial asymmetry and 
oq is the Gaussian depolarization rate). The exponen- 
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FIG. 1: (Color online) The representative ZF and TF-/xSR 
time-spectra measured above and below T c (left panels), and 
the dependence of the superconducting part of square root 
of the second moment a BC oc A~ ft 2 oc p s on temperature and 
magnetic field (right panels) for UD23K (a), OP35K (b), and 
OD29K (c). 



tial character of the muon polarization decay in UD23K, 
described as A ZF (t) = Aoexp(-At), is an indication of a 
weak magnetism which is probably caused by the closed 
proximity to the static magnetic order fl6j . 

The TF-/xSR data were analyzed by using a two- 
component Gaussian fit of the /iSR time-spectra allow- 
ing to describe the asymmetric local magnetic field dis- 
tribution P{B) in a superconductor in the vortex state: 

A™(t) = TLi M exp(-<7?i 2 /2) cos^i + if) 0, E3- 
|l9j. Here Ai, <7j, and Bi are the initial asymmetry, re- 
laxation rate, and mean field of the i— th component, 
7 M = 27r x 135.5342 MHz/T is the muon gyromagnetic 
ratio, and ip is the initial phase of the muon-spin ensem- 
ble. The weak magnetism detected in ZF experiments for 
UD23K was taken into account by multiplying the fitting 
function by exp(— At). The superconducting part of the 
square root of the second moment a sc oc A~ b 2 oc p s (X a b is 
the in-plane magnetic penetration depth) was further ob- 
tained by subtracting the normal state nuclear moment 



contribution (<r n 
P(B) (a 2 ) as cr s , 



from the measured second moment of 
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The magnetic field dependence of a sc at T = 1.6 K for 
the samples studied is shown in the corresponding insets 
of Figs. Q] a, b, and c. The decrease of a sc for OP35K 
and OD29K is a consequence of both, the nonlinear and 
the nonlocal response of the superconductor containing 
nodes in the energy gap to the increasing magnetic field 
[ilj . The solid lines for OD29K and OP35K correspond 
to fits of the relation a SC {H) / a SC (H = 0) = 1 — K- \/H to 
<j sc (H) which takes into account the nonlinear correction 
to p s for a superconductor with a d— wave energy gap 
[20l | . The analysis of cr sc (H) for UD23K reveals, however, 
that only at very low fields (< 30 mT) a sc follows the 
tendency observed for OP doped and OD Bi2201 samples. 
For higher fields a sc increases with increasing H. Such 
behavior is generally associated with the field induced 
magnetism and is often observed in various underdoped 
cuprate HTS (see e.g. Refs. ED, 



The observation of the field induced magnetism in 
UD23K for fields exceeding 30 mT is quite unexpected, 
especially considering the fact that in their recent study 
Russo et al. [l6[ do not detect any kind of field induced 
effects in Bi2201 with T c ~ 27 K up to p Q H ~ 5 T. We 
may suggest, therefore, that the partial substitution of 
Bi by Pb, as made in our samples, leads to enhancement 
of coupling between Cu02 layers, thus causing Bi2201 
to be more 3-dimensional. This could affect both, the 
superconducting and the magnetic properties. As for su- 
perconductivity, our previous studies point to a substan- 
tial reduction of the anisotropy coefficient 7A = X c /X a b 
(A c is the out-of plane component of the magnetic pene- 
tration depth) as well as to the shift of the vortex-lattice 
melting transition much closer to T c 14|. The magnetic 



properties could be also affected due to increase of the 
interlayer magnetic exchange coupling J', which for the 
layered materials like cuprate HTS is the primary quan- 
tity determining the Neel temperature Tn oc J' £| D (£20 
is the magnetic correlation lengths within the layer). 

The temperature dependence of X a b was obtained from 
the measured cr sc (H = const, T)'s shown in Fig. [T] 
by following the procedure described in Ref. [l5|. It in- 
cludes, first, the reconstruction of the effective penetra- 
tion depth X e f f (H, T) from a sc (H, T) by using the rela- 
ti on a sc {H,T ) = 4.83 • 10 4 [1 - H/H c2 (T)][l + 1.21(1 - 
y/H/H c2 {T)f] X~f } [11]. Here H c2 is the upper critical 
field with the zero-temperature values (J>oH c2 {0) — 60 T, 
~ 50 T and ~ 45 T for UD23K, OP35K and OD29K, 
respectively [24| and with the T dependence follow- 
ing the Werthamer-Helfand-Hohenberg prediction pol ]. 
As a next step, A a & was reconstructed by decomposing 
X e f f (H, T) into the field and the temperature dependent 
components as X eS (H/ H c2 ,T) = C{H/H c2 ) X ab {T). 

The dependence of X~ b 2 oc p s on temperature is shown 
in Fig. H For OP35K and OD29K A~ b 2 (T) was recon- 
structed from a sc (H,T)'s measured at hqH = 0.04, 0.1, 
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FIG. 2: (Color online) (a) The pSR and ARPES data for 
UD23K sample. Left panel: A~ 6 2 (T) oc p s (T). Lines are the 
simulated A~ b 2 (T) curves obtained by assuming that the full 
spectral gap measured by means of ARPES leads to the pair- 
ing and QW(0) = 1 (green line) and for the gap of d— wave 
symmetry with QW(</>) measured by ARPES (red line). Top 
right panel: angular dependence of the quasiparticle weight 
function QW(<^>) [jlj . <f> — and <j> = 7r/4 correspond to the 
antinodal and the nodal directions, respectively. Bottom right 
panel: angular dependence of the spectral gap above (41 K, 
red points) and below (11 K, blue points) the transition tem- 
perature T c — 23 K [Tjj]. The solid line is the cos 20 fit of 
T — 11 K data in the nodal region, (b) and (c) - the same as 
in (a), but for OP35K and OD29K, respectively. 



0.2, 0.4 and 0.64 T (see Ref. QJ), and p Q H = 0.04 and 
0.2 T, respectively. The field induced magnetism does not 
allow us to make a similar reconstruction for UD23K. We 
assumed, therefore, X a b(T) ~ A e //(0.01 T, T). As shown 
in Ref. [l5| the relation A e //(ff, T) ~ A a t,(T) is still valid 
in the limit of low magnetic fields H < 10~ 3 -ff C 2(0). 

The resulting A~ fc 2 (T)'s were a naly zed within the same 
scheme as described in Refs. Il4 1 15 : 
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FIG. 3: (Color online) Temperature dependence of the su- 
perconducting gap for UD23K (a), OP35K (b), and OD29K 

(c) . The open and closed symbols correspond to QW(4>) = 1 
and QW(0) measured by means of ARPES [lit- Tne s °l id 
lines are BCS fits with parameters summarized in Table [T] 

(d) Dependence of A~ t 2 (T = 0) and the coherence peak inten- 
sity Wcp integrated over the whole Fermi surface on T c . (e) 
Angular dependence of Wcp measured by ARPES [TlT |. 



the angle along the Fermi surface (see the top right panels 
in Figs.[2]a, b, and c), A(T, </>) denotes the superconduct- 
ing (pairing) gap depending on T and 4>, QW(0) accounts 
for the relative weight of the quasiparticles condensed 
into the Cooper pairs, and Sqw = Jq^ 4 QW(</>)d</>. Note 
that we do not consider here the presence of the second 
s— wave gap 3, 2(| 27 1 since the d— wave contribution 
to A 



lb (T) in Bi2201 seems to be predominant. 
As a first step we checked if the full spectral gap 
measured in ARPES experiments (lower right panels in 
Figs. [5] a, b, and c) could be a pairing gap. The solid 
green lines on the left panels of Figs. [2] a, b, and c cor- 
respond to the theoretical curves obtained by assuming 
that the states over the whole Fermi surface are equally 
available for condensation [QW(</>) = 1], and that the 
gap in the antinodal region is T-indcpcndent and it fol- 
lows the BCS prediction close to the nodes [l(j. The 
poor agreement between the experiment and the theory 
suggests that the gap near the antinodes could not be 
related to the pairing. In particular, the step-like jump 
of p s at T — T c is due to the fact that the antinodal gap 
is not closed when the temperature passes through T c . 

We suggest, therefore, that for all levels of doping the 
superconducting (pairing) gap has a d— wave symmetry: 
A(T,<j>) = A(T) cos 20. This is indeed true for OD29K 
(Fig. [2] and Ref. [ll|), as well as for various OD and OP 
hole-doped cuprate HTS studied by means of ARPES 
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(see e.g. Ref. l28l and references therein). The predomi- 
nantly d— wave symmetry of the order parameter in UD 
HTS was also confirmed in tricrystal experiments p9| . 
This, together with known QW((f>) and A~ h 2 (T) allows 
the use of Eq. (fl]) to reconstruct A(T). The results of 
such reconstruction for two possible scenarios are shown 
in Fig. [3j the first with QW(j£) obtained by means of 
ARPES (see Fig. El and Ref. [3), and the second with 
QW(0) = 1, indicated by the closed and open symbols 
respectively. Fits of the BCS model to A(T) are repre- 
sented by solid lines. Both sets of A(T) data follow the 
BCS temperature dependence in agreement with the re- 
sults of ARPES studies for the spectral gap in the nodal 
region (lo| . 



TABLE I: Zero-temperature values of the superconducting 
gap A(0) and the ratio 2A(0)/fcsT c as obtained from the fit 
of the BCS model to A(T) represented in Figs. [3] a, b, and 
c. The "d-wave+ QW((j>) = 1" and "d-wave+QW(^) " refer 
to the case of QW(0) = 1 and QW(0) measured by means of 
ARPES [ll[, respectively (see text for details). 



d-wavc+ QW(4>) = 1 

A(0) m 
(meV) 


d-wavc+QW(c» 
A(0) 2A(0) 

k B T c 

(meV) 


UD23K 4.6(2) 4.6(2) 
OP35K 9.3(2) 6.2(1) 
OD29K 8.6(2) 6.9(2) 


8.4(3) 8.5(3) 
13.6(3) 9.0(2) 
10.3(3) 8.2(2) 



In order to distinguish between the two above men- 
tioned scenarios we note that: (i) Accounting for the 
quasiparticle weight as measured by ARPES causes a sys- 
tematic shift of A(T) to higher values. This leads to bet- 
ter agreement with the gap obtained from the cos 2c/> fit to 
the ARPES data in the nodal region ll| . (ii) The ratio 
2A/fcsT c increases with doping from ~ 4.6 for UD23K to 
6.9 for OD29K in a case when QW(</>) = 1 while it stays 
almost constant (~ 8.5) for QW(<p) obtained by means 
of ARPES, see Table Q] Note that the independence of 
2A/ksT c ratio on doping is well confirmed experimen- 
tally for various HTS families (see e.g. Refs. 30, 31 and 
references therein), (iii) The intensity of the coherence 
peak Wcp (Fig.[3]e) integrated over the whole Fermi sur- 
face scales with the zero-temperature superfluid density 
(Fig. |3]d), thus pointing to the direct relation of Wcp to 
the " local " superfluid density 0, 0, EH ■ By approaching 
the antinodal point, Wcp decreases thus requiring the 
corresponding decrease in the local density of the super- 
carriers. The strongest effect (decrease of Wcp down to 
0) is observed for UD23K while for OD29K the weight 
of CP in the nodal region still remains substantial. All 
these arguments taken together support the scenario ac- 
cording to which the gap (pseudogap) near the antinodes 
makes a part of the states unavailable for the supercon- 



ducting condensation and leads, therefore, to a reduced 
density of the supcrcarriers in the antinodal region. 

To summarize, the temperature dependence of the su- 
perfluid density p s was studied in underdoped (T c ~ 
23 K), optimally doped (T c ~ 35 K) and overdoped 
(T c ~ 29 K) single-crystalline Bi2201 samples by means 
of muon-spin rotation. By comparing the measured 
p s (T) with that calculated theoretically based on the re- 
sults of ARPES [lH we found that the superconducting 
gap in Bi2201 at all levels of doping has d— wave sym- 
metry and that A(T) follows reasonably well the BCS 
prediction. It was also shown that p s (T) is inconsistent 
with the case when the carriers over the whole Fermi sur- 
face arc equally available for condensation thus suggest- 
ing that some parts of the Fermi surface do not develop 
the superconducting coherence. 
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